Graphene quantum dots formed by a spatial modulation of the Dirac gap 
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An electrostatic quantum dot cannot be formed in monolayer graphene, because of the Klein 
tunnelling. However, a dot can be formed with the help of a uniform magnetic field. As shown 
here, a spatial modulation of the Dirac gap leads to confined states with discrete energy levels, thus 
defining a dot, without applying external electric and magnetic fields. Gap-induced dot states can 
coexist and couple with states introduced by an electrostatic potential. This property allows the 
region in which the resulting states are localized to be tuned with the potential. 



The relativistic character of electrons in graphene 
has attracted considerable attention [l|, This work 
presents a way of forming a graphene quantum dot as a 
consequence of a spatial modulation of the Dirac gap. 

In an ideal graphene sheet, the band structure has no 
energy gap and the Dirac electrons are massless. The 
quantum states of a graphene dot, formed by an ex- 
ternal electrostatic potential, are deconfined due to the 
Klein tunnelling that is inherent to massless particles Q . 
Therefore, an electrostatic dot cannot confine electrons 
in graphene since these can tunnel through any potential 
barrier. However, a uniform magnetic field suppresses 
the Klein tunnelling leading to confined states 0, Q . 

It is experimentally possible to engineer an energy gap 
in graphene's band structure, referred to as a Dirac 
with a value ranging from a few to hundreds of meV [ 
The gap leads to electrons with mass and thus an electro- 
static potential results in confined states as in common 
semiconductors [11 . Most importantly, as shown here, a 
spatially-modulated gap induces confined states regard- 
less of the application of external fields. This can be 
achieved provided the gap has a local minimum in which 
the states become localized, thus defining a quantum dot. 
The application of a potential, generated by a gate elec- 
trode, couples the gap- and potential-induced states. The 
coupling strength is tunable with the potential, and it 
determines the region in which the resulting states are 
localized. 

A spatially- modulated Dirac gap has been reported Q . 
Possible ways of creating the required gap modulation 
that forms the quantum dot include substrate engi- 
neering and the application of strain to the graphene 
sheet [i|,l3. 

The physics of a graphene dot, for energies near the 
Dirac points, is described by the Hamiltonian Q 



H = vpcr ■ (p + eA) + VI + rAcr^ 



(1) 



where the Fermi velocity vp=^ /h, with 7=646 nieV nni, 
is assumed position-independent. cr={ax,(Jy), (Jz are the 
2x2 Pauli operators acting on the two carbon sublattices, 
p=—ifiV——ih{dx, dy) is the 2D momentum operator, A 
is the vector potential that generates the magnetic field 
B=VxA, V is the electrostatic potential and I is the 
unit matrix. The last term in Eq. ([T]), referred to as mass 
term, gives rise to an energy gap 2A in the spectrum of 



graphene, where r=l (t=— 1) corresponds to the K [K') 
valley. 

For the dot model, V and A are chosen cylin- 
drically symmetric and the magnetic field is uniform 
and perpendicular to the graphene sheet, B=i?z, so 
A=(0,Ae,0), with Ae=Br/2. The Dirac equation 
H^=E^ , can be written in cylindrical coordinates with 
^'=r~^/^{/i(r) exp[i(m — 1)6*], i/2(r) exp(«r7i6')}, where 
m=0, ±1... is the angular momentum quantum number. 
The radial functions /i and J2 satisfy 



(2) 
(3) 



0, 



U--i-]h + {V-E- TA)f2 



with [/=7(2to- l)/2r + 7eBr/2;i. Equations © and dH) 
are satisfied for both confined and deconfined states. The 
former have an exponential tail asymptotically, e.g., in 
the limit of large radial distance r (r— >-oo), whereas the 
latter have an oscillatory tail. If for large r V and A 
are constant or have a power-law dependence, then the 
confined-deconfined character of a state is determined by 
the asymptotic sign of 
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(4) 



A state with energy E is confined only if q is asymp- 
totically negative Otherwise the state is deconfined. 
This criterion indicates that confined states can be in- 
duced even for B=0 and V=0 everywhere, provided that 
E^—A^<0 asymptotically. This inequality cannot be sat- 
isfied when A is spatially-independent because all the en- 
ergies satisfy |i?|>A. But the inequality can be satisfied 
when A is spatially-dependent with an asymptotic value 
larger than that for small r. This happens, for example, 
when A is zero within a disc area, and nonzero outside 
that area, A=So. Then a number of discrete energy lev- 
els satisfy |i?|<(5o and correspond to confined states with 
a large amplitude within the disc area. These states can 
be regarded as dot states. 

In the presence of an electrostatic potential V and 
B=0, Eq. Q shows that if V and A are unequal and 
rise asymptotically, then confined states occur only if 
{V — A) <0 so that q<0. In this case confinement is 
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FIG. 1: Energy levels versus the asymptotic value of the mass 
term So, that generates a Dirac gap 25o. 



energy independent. However, this work focuses on the 
case where both V and A are constant asymptotically, 
which is the most common experimental regime. Then 
confinement occurs if {V — E)'^ — A^<0 and thus it is 
energy-dependent . 

The properties of the gap-induced dot are analyzed by 
solving numerically the two coupled equations [Eqs. ([2]) 
and ([3])] using a discretisation scheme, which satisfies 
time-reversal symmetry: E(rn, B)=E{1 — m,—B) for 
A=0 and E(m, B, t)=E{1 - m, -B, -r) for A^^O. The 
spatially-dependent mass term is modelled by A=0 for 
r<R and A(r)=-(5o/cosh[(r--i?)/rf]2-f(5o for r>R so that 
asymptotically A«(5o- This choice is not of particular im- 
portance; either a smooth or sharp modulation of A re- 
sults in confined states. For brevity, all the results shown 
are for i?=250 nm, c?=150 nm, t=1 and m=5. 

Figure [T] shows the energy levels as a function of the 
asymptotic value of the mass term 5q for B=Q and V=Q. 
Deconfined states correspond to the (quasi) continuum of 
levels indicated by the black area in Fig. [T] [9] . Confined 
states correspond to the discrete levels emerging through 
the continuum via anticrossing points (see below), form- 
ing two distinct ladders of energy separated by a gap. 
The spacing of the discrete levels and the gap increase 
with ^0- Further, the confinement becomes stronger with 
(5o, and therefore an increase in. 5n leads to an increase in 
the number of discrete levels However, for a fixed 

(5o, not all m values give confined states, since the angu- 
lar momentum term U , for B=0, tends to delocalize the 
states for large m. For this reason there are no confined 
states for 6o<10 meV in Fig.[T] The physics is different if 
B^O and V is constant since from Eq. (jj]) q<0 asymptot- 
ically, leading to confined states independent of energy. 
One such case is illustrated in Fig. [T] for 5=0.1 T. The 
continuum of levels has been replaced by a discrete set for 
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FIG. 2: (a) Energy levels versus the potential depth Vb for a 
constant Dirac gap 2A=50 meV. (b) As (a) but for a spatially- 
modulated Dirac gap with an asymptotic value of 25o=50 
meV. The states of the energies marked by circles are shown 
in Fig. O Enlarged views of (a) and (b) are shown in (c) and 
(d). 



all do, reflecting the disappearance of deconfined states. 
This is valid for all m values and r=— 1. 

The effect of an electrostatic potential on the gap- 
induced dot is now investigated. The potential that is 
generated in the graphene sheet by gate electrodes can 
be calculated within the Thomas-Fermi model A 
slowly varying quantum well potential is approximated 
by y(r)=— Vo exp(— t^/Zq). The quantum well depth is 
Vq, the width is Iq, and asymptotically V~0. 

Figure [2] shows the energy level diagram as a function 
of Vb for ^0=180 nm. Confined and deconfined states are 
identified as in Fig. [TJ For a constant gap and small Vq, 
the angular momentum delocalizes the states; therefore 
confined states are formed after a critical value of Vq . The 
general trend is that with increasing Vb, the number of 
discrete levels increases while the lowest levels merge into 
the continuum. When this happens, the corresponding 
states undergo a transition from confined to deconfined, 
which is refiected in the energy diagram by the appear- 
ance of anticrossing points [Fig. [5](c)]. These also appear 
when the states undergo the opposite transition for en- 
ergies near Sq. For a spatially-dependent gap, there exist 
discrete levels even for Vb=0, because of the gap-induced 
confinement. Unlike the constant-gap system, as Vb in- 
creases, anticrossing points are formed between discrete 
levels [Fig. [2l^d)], reflecting a coupling between confined 
states due to the potential and the spatially-dependent 
gap. This coupling is strong for the gap-induced states of 
the upper ladder of energy and therefore the correspond- 
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FIG. 3: Quantum states for different potential depths Vo, 
from top to bottom: Vb=0, 30, 75, 120 meV. The Dirac gap 
is spatially modulated with an asymptotic value of 25o=50 
meV. Left (right) panels show states with energies marked by 
• (o) in Fig. [lib). The vertical axes of the insets range from 
to 2x10"^. 

ing anticrossing points are not well-defined. In contrast, 
states of the lower ladder couple weakly to the potential. 

Figure |3] illustrates the effect of the potential on dot 
states, for the energies shown in Fig. [Hb). For Vo=0, 
the states are confined owing to the gap modulation. 
Consider now the states shown in the left panels. As 
Vb increases, the gap-induced state with positive energy 
couples to the potential, e.g. for Vb=30 meV, and with 
increasing Vb it becomes localized in a region defined by 
the potential. This state then couples to gap-induced 
states of the lower ladder, e.g. for Vb=75 meV, while 
its energy decreases. For Vb=120 meV the state is de- 
confined, with an oscillatory tail, and its energy lies in 
the continuum. A coupling between gap- and potential- 
induced states occurs also for the states shown in the 



right panels, e.g. for T4)=75, 120 meV. As Vb increases 
the state with the maximum energy in the lower ladder 
(for Vb=0), couples with excited potential states. This 
trend is consistent with the series of anticrossing points 
in the Eiy^) plot. 

In Fig. [3] the coupling between potential- and gap- 
induced states of the lower ladder is weak (Vb=75 meV); 
therefore the states peak in the region defined either 
by the potential profile (left) or the gap modulation 
(right) These two regions have a small overlap when 
R>1{) and Vb is large. Strong coupling can be induced for 
small m; for instance, for m—X the states peak strongly 
in both regions. The m-dependence of the coupling can 
be explained within a semiclassical approach The 
relative maximum amplitude of the two components, as 
can be derived from Eqs. ^ and satisfies |/i|~|/2| 
within a A=0 region, whereas |/i|>|/2| (|/i|<|/2|) for 
energies in the upper (lower) ladder. The latter behavior 
is more pronounced when A is large and constant with 
Vb^^O. Then one of the components becomes vanishingly 
small depending on choice of energy and valley (t=±1). 

In order to probe the states of the gap-induced dot 
the Fermi level has to be adjusted near the middle of the 
gap, where only confined states with small values of m lie, 
and hence the resultant density of states is low. For the 
same reason the electrostatic potential has to be small. 
Then, it should be experimentally possible to resolve the 
quantum states using similar measurements as in GaAs 
quantum dots. 

In summary, a graphene dot can be formed as a result 
of a spatial modulation of the Dirac gap, without apply- 
ing external fields. An electrostatic potential allows gap- 
and potential-induced states to coexist and become cou- 
pled as the potential increases. The coupling strength 
determines the region in which the states are localized. 
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